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Abstract 

This progress report consists of two separate reports. The first one describes 
our work on the use of variable gain amplifiers to increase the receiver dynamic 
range of space borne laser altimeters such as NASA’s Geoscience Laser Altimeter 
Systems (GLAS). The requirement of the receiver dynamic range was first 
calculated. A breadboard variable gain amplifier circuit was made and the 
performance was fully characterized. The circuit will also be tested in flight on 
board the Shuttle Laser Altimeter (SLA-02) next year. 

The second report describe our research on the master clock oscillator 
frequency calibration for space borne laser altimeter systems using GPS 
receivers. 
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1. Introduction 

The received signals of a space borne laser altimeter varies over a 
broad range due to variation in the reflectivity and slope of the 
ground targets. The target reflectivity affects the received signal 
pulse energy, the slope affects the pulsewidth, and both affect the 
pulse amplitude. The photodetector and the subsequent amplifiers 
have to be able to accommodate the dynamic range of the received 
pulse amplitude, with sufficient gain to overcome the preamplifier 
noise but without saturation at anytime. Since the signal processing 
circuit such as an A/D converter usually has a limited input dynamic 
range, a variable gain amplifier (VGA) is necessary to achieve the 
optimal overall receiver performance. 

We are currently developing a photodetector assembly for the 
Geoscience Laser Altimeter System (GLAS), which will include a 
wideband and wide control range VGA. The up coming Shuttle Laser 
Altimeter II (SLA-II) also requires a VGA in the photodetector 
assembly to improve the receiver dynamic range. We have built an 
engineering model VGA circuit card for SLA-II to substitute the 
first post amplifier card inside the detector assembly. Our test 
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results show that the circuit meets all the requirements for SLA-II. 
Its performance also exceeds the VGA requirement for GLAS. 

The rest of this report describes the details of the VGA requirement 
for GLAS and SLA-II and a proposed GLAS detector assembly circuit 
design which includes a photodetector preamplifier module, a VGA, a 
low power amplifier, and a power amplifier. Lastly, the design and 
the test results of the SLA-II VGA circuit card will be presented. 

2. Signal Dynamic Range Requirement for GLAS 

The expected signal level at the GLAS altimeter receiver is 1 .45 fJ 
(7760 photons) per pulse incident to the photodetector, assuming 80 
cm diameter telescope, 25 cm diameter center obscuration, 100 mJ 
laser pulse energy, 20% target reflectivity, 705 km altitude orbit, 
50% optics transmission, and 50% round trip atmosphere 
transmission 1 . The peak received signal power is given as the 
received pulse energy divided by the pulse width. The most likely 
received pulse width is 30 ns FWHM, which corresponds to a 3 degree 
slope ground target and 70 m laser footprint. The received peak 
power under this pulse width is, therefore, 1.45fJ/30ns = 48.3 nW. 

The peak signal power can increase by 24 times (27.6dB) when GLAS 
is over flat ice sheet with 80% reflectivity and 5 ns received pulse 
width. Note it is the signal peak power that we have to accommodate 
without saturation and to give an output signal which is within the 
optimal range for the subsequent signal processing circuit. 


1 see the attached link margin analysis. 
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The lowest signal peak power to be considered should be the same as 
the minimum detectable signal, which is about 130 photons/pulse 
for 30 ns FWHM at nighttime. The minimum detectable signal level 
also increases with the received pulse width, but we consider only 
30 ns here as the typical case. The overall input signal range to be 
considered is, therefore, 7760/130 times 24, or a factor of 1400 
(63dB) for nighttime. As a comparison, the range for daytime is a 
factor of 470 (53dB). 

To achieve the optimal receiver performance, the VGA gain 
adjustment range plus the optimal input signal range for the 
subsequent signal processing circuitry such as a waveform digitizer 
should be greater than or equal to the maximum input signal range at 
nighttime. For a 8 bit waveform digitizer, the optimal input signal 
lies between 1/16 to 1/2 the full scale (4-7 bits out of the full 8 
bits) based on our analysis and the experimental results 2 , which 
gives an optimal input range of 8, or 18 dB. One can further narrow 
the required VGA range by another factor of 2 (6 dB) because the 
signal level into the waveform digitizer can be lower than the 
optimal level for weak signal when other noise sources dominate 
that 3 bits A/D resolution is sufficient 3 . Therefore, the VGA range 
can be as low as 39 dB. 


2 X. Sun, “Design and Performance Measurement of the GLAS Breadboard Altimeter 
Receiver," Interim Progress Report on NASA Grant NAG5-2232, "Receiver Design, 
Performance Analysis, and Evaluation for Space-Borne Laser Altimeters and Space-to- 
Space Laser Ranging Systems," for the period of Oct. 16, 1995 to April 14, 1996, Dept. 
ECE, The Johns Hopkins University, Baltimore, MD 21218. 

3 Same as above. 
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Another 3 dB reduction in the VGA range could be obtained if one 
designed the photodetector preamplifier such that the impulse 
response matches the received optical signal pulse that the output 
pulse width is V2 that of the input. However, this is difficult to 
achieve in practice, because the primary goal of the preamplifier is 
low noise and wide band, preferably not to be complicated by the 
impulse response pulse shape. 

The electrical bandwidth of the VGA and the subsequent amplifiers 
have to be wide enough to pass the shortest possible received laser 
pulses, 5 ns FWHM. The APD preamplifier module has a 3 dB 
bandwidth of about 150 MHz with a RC filter type roll-off. The 
subsequent amplifiers should all have similar or wider bandwidth. 

The gain of the photodetector preamplifier hybrid module for GLAS 
is about 200 kV/W. The output dynamic range is greater than 60 dB 
according our previous test results. The output pulse amplitude in 
response to the 1.45fJ /30 ns pulse is roughly 10 mV. The optimal 
input signal range for a waveform digitizer is typically about 
0.5 volt (mid scale). Usually the signal output from the detector 
assembly is split between the waveform digitizer and other 
branches of the receiver such as the time interval unit (TIU). The 
desired signal amplitude from the detector assembly should be about 
1 volts. Consequently, the design goal for the detector post 
amplifiers including the VGA should be 40 dB nominal gain, variable 
by at least +12 to -28 dB (40dB range), with a roughly 10 mV 
nominal input and 1 volt output pulse amplitude (13dBm into 50Q). 
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Usually a VGA and several linear amplifiers have to be used in 
cascade to achieve this goal. 

3. VGA Requirement for SLA-II 

SLA-II is similar to GLAS except for the lower orbit altitude, 310 
km, smaller telescope diameter, 40 cm, lower transmitted laser 
pulse energy, 40 mJ/pulse, and wider transmitted laser pulse width, 
10 ns. The expected received signal level for 0.20 reflectivity 
target is 0.52 times that of GLAS, or 0.75 fJ/pulse (4040 photons). 
The peak optical signal power for 30 ns FWHM pulses is 25 nW. The 
maximum input signal is 12 (21.6dB) times higher when it is over 
flat ice or snow ground target (0.8 reflectivity and 10 ns pulse 
width). The minimum detectable signal is roughly the same as 
SLA-II (130 photons/30ns pulse). The input signal dynamic range is, 
therefore, 370 or 51 dB. The optimal input signal range for the 
waveform digitizer and TIU should be the same as GLAS, 18+6=24 dB. 
Therefore, the VGA range required is about 27 dB. 

The 3dB bandwidth required for the SLA detector is greater than 
75 MHz to pass 10 ns FWHM pulses. 

The control signal needs to be 0-5 volts since a D/A control signal 
from the on board computer is readily available in the current 
system. The DC power supplies are preferred to be ±12 V. 

The APD preamplifier hybrid module has a gain of 775 kV/W. The 
output pulse amplitude at 30 ns FWHM is 19.4 mV. The optimal 
output signal amplitude from the second stage should be about 
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2 volts, which leads to a total of 40 dB gain of the entire detector 
assembly. At present, SLA detector assembly has two post amplifier 
stages, one has a voltage gain of 24.2 (27.7dB) and the other 10 
(20dB). We can just replace the first amplifier stage with a VGA 
circuit card and leave the second stage unchanged. Consequently, we 
need a VGA circuit card of 20 dB nominal gain and variable by +5 to 
-22 dB with a nominal input signal of about 20 mV. 

4. Proposed GLAS Detector Post Amplifier Circuit 

One candidate VGA for both GLAS and SLA is the wideband AGC 
amplifier, Model SL6140 by GEC Plessey Semiconductors 4 . It has a 
gain from at least -32 to 15 dB, a 53 dB range with 50 Q load (70 dB 
specified for 1KQ load), 400 MHz bandwidth, and a maximum output 
amplitude of 600 mV into 50 Q load (8.6dBm). 

One candidate for the last stage is wideband low distortion power 
amplifier, CLC561 by Comlinear Corp. It has a 3 dB bandwidth of 150 
MHz at +24dBm output level (10V peak to peak). The gain can be set 
to as high as 26 dB without significant reduction in the bandwidth. 

To achieve the required gain and variation range, 40 dB plus 12 
minus 28 dB, one can set the nominal VGA gain to 0 dB and use 
another 14 dB low power amplifier before or after the VGA. This 
type of linear amplifiers are relatively easy to find. One candidate is 
the GEC Plessey SL560 or its equivalence, which has a bandwidth of 
300 MHz and a gain of 14 dB 5 . 


4 see the attached data sheets. 

5 See the attached data sheets. 



5. Engineering Model VGA Circuit Card for SLA-II and the 
Test Results 

We have built an engineering model VGA circuit card for SLA-II, to 
replace the post amplifier #1 card inside the detector assembly. The 
circuit consists of a GEC Plessey SL6140 VGA amplifier followed by 
an SL560 linear amplifier. Figure 1 shows the circuit schematic. The 
circuit is actually part of the amplifier chain which we proposed for 
the GLAS detector assembly in the previous section. The circuit was 
fully tested and well met all of our design goals. 

The entire circuit was operated from a single +12V power supply. 
The SL560 amplifier required a lower supply voltage and a 150 Q 
resister was used to drop the 12 volt supply to the desired 6-9 volts 
range. The current drawn by this amplifier was 30 mA. The VGA 
stage drew 18 mA. The total power consumption of the card was 
580 mW. 

Figure 2 shows a block diagram of test setup. The input signal 
consists of 20 ns wide rectangular pulses. The oscilloscope 
monitored simultaneously both the output and the input (times 10 
due to the 20 dB pad) and measured the pulse parameters. The VGA 
control voltage was generated by a programmable power supply. 

Figure 3 gives a spreadsheet of the test conditions and the results. 
Figure 4 plots the voltage gain and the relative propagation delay vs. 
control voltage. We have achieved a gain control range from -23 to 
30 dB, well exceeded the range required by SLA-II (-22 to 5 dB). The 
propagation delay varied by no more than 0.5 ns over the entire gain 
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range. The maximum output pulse amplitude was found to be about 
600 mV. 

Figure 5 shows the input and the output pulse at close to the 
maximum gain. It shows the circuit has a rise and fall times within 
2.4 ns. There was little pulse overshoot and undershoot. Figures 6 
through 9 show the measured input and output pulses at maximum, 
10, 1, and 0.1 voltage gain settings. At very low gain, the output has 
some overshoot and undershoot, as shown in Figure 9. 

Figure 10 shows the input and the output pulse shapes at individual 
amplifier stages. The VGA stage output looked noisy because of the 
oscilloscope probe. The gain of the linear amplifier (called "buffer") 
after the VGA is shown to be 14.5 dB. The gain of the VGA stage is, 
therefore, -37.5 to 15.5 dB. 

Figure 11 shows the slight pulse undershoot due to AC coupling for 
relatively wide pulses, 1 [is in this case. The AC coupling capacitors 
were 0.47 |iF. Since the input and output impedance were 50 ft , the 3 
dB lower cutoff frequency was about 6.8 kHz, which was the same as 
the rest of the detector assembly. The lower cutoff frequency can be 
reduced by increasing the capacitance or using two or three of those 
capacitors in parallel. 

We also measured the VGA response time to the control signal by 
feeding it with a 2 MHz pulse train and applying a step signal to the 
control terminal. The bypass capacitor, C3, was temporally removed 
during this measurement. Figures 12 and 13 show the measurement 
results. The response time of the VGA was less than 0.5 jis without 


8 



any ringing. Therefore, we are capable of changing the detector 
assembly gain on a pulse to pulse base for both SLA-II and GLAS. 
Currently, the VGA response time was set to about 1 ms with the 
use of a 0.1 |iF bypass capacitor. 

Figures 14,15, and 16 show the S-parameter measured with our 
HP8753C/85047 network analyzer from 300 kHz to 500 MHz. The 
3 dB bandwidth of the circuit is shown to be about 300 MHz and the 
phase response was very close to linear. The input voltage standing 
wave ratio (VSWR) was less than 1.25:1, mostly because of the 50 ft 
resistor at the input. The output stage inside SL560 was an emitter 
follower, which had a very low source impedance and the measured 
S 22 was not shown here. 

Lastly, Figure 17 shows the measured quiescent output noise 
spectrum measured by the HP3585A and HP8561B spectrum analyzer. 
The average spectral density was -128 dBm/Hz, which corresponded 
to 1.54 mV rms on a 50 ft load over the 300 MHz bandwidth, which 
was about what we saw on the oscilloscope during the test. The 
highest noise spikes within the 300 MHz bandwidth were -65 dBm at 
about 20 MHz, which corresponded to an insignificant 0.1.3 mV rms 
in time domain. Note the circuit card was not well shielded during 
this test and the spikes might have come from RF pick-up from the 
environment. 
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GLAS Receiver Pulse Width Calculations 
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GLAS Receiver Pulse Width Calculations 
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N := 100 At = .1 

.. N 

t = i- At At 

1 1 


1 = 0..N 


Y i '* y (*i) F i = f (^) Y max = ““(Y) f max = ^(F) 



RMS pulse width 


4 


y (t.) 2 - y. 

^ 1 \ ■/ * 


VV 1 =2.104 


FWHM pulse width 

k =0.100 Zl k = F n k tl k =t N _ k j = 0.. (N - 10) Z2j = 


ENFWHM = 2 linterp Z l , z 1 , 


OUTFWHM -- 2 


f \ 
1 max 1 

2 
Y 


linterp! Z2,x t - 


max 




1NFWHM =5 


OUTFWHM =5 


"n-J T j = Sj - j 




N = 150 At =3 i - 0..N 

t. -- i- At - -~*At- Y; =y(t.,a) Y ^ = max(Y) 



RMS putsewidth 
FWHM pulse width 
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GLAS Altimeter Link Margin Analysis 

Xiaoli Sun, The Johns Hopkins University 

System parameters 


E t := 0.0740 

Laser transmitter pulse energy (J) 

R:= 598 10 3 

Range to the target (m) 

r diff := °- 60 

Target diffusion coefficient 

♦tel :=0 - 90 

Receiver Telescope Diameter (m) 

♦ obs :=0 - 25 

Telescope center obscuration (m) 

^ rcvr ;= 0 5 

Receiver optics transmission 

9 POV := 0-^0' 10 

Receiver FOV, full angle (rad) 

Ak := 0.60- 10' 3 

Receiver optical bandwidth (um) 

0 atmo := 0 5 

Atmosphere Transmission Coefficient (one way) 

1 solar := 

Solar irradiance at the target (W/m A 2/um) 

Other system parameter 

values : 

X:= 15.242- 10' 9 

Received laser pulse width (receiver integration time) 

^ APD := 035 

APD quantum efficiency 

G := 150 

Average APD gain 

k eff := 001 

APD ionization coefficient ratio 

I b := 50- 10" 12 

APD bulk leakage current (A) 

I s := 15- 10' 9 

APD surface leakage current (A) 

N a,.ip : = t 6Kf T 

Preamp spectral noise current density (A A 2/Hz) 

2-jct 

Receiver noise bandwidth (Hz) BW n = 2.30 1 4* 10 7 

O 

© 

rv» 

ii 

o 

Range gate (km) 
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RECEIVED SIGNAL AND NOISE 

Solar background noise photons per second 


P b := 1 solar'^ ^ rcvr' 71 ’ 
^ APD' P b 


Iq fov^ 

2 I 

r diff 

2 f 

K \ 


1.17 1.6 10 




19 


n b := ^b' T 


P b = 1.7119-10 ' n b =48.79 

Received signal pulse energy 

_ _ r diff | ♦ tel " ♦ obs I 2 

E r : = E t 1 *' ; h fcvr’^t 

K \ 4R 

Signal photons per received pulse 
E r 

N s :=- 


1 rcvr ** atmo 


1.17* 1.6* 10 


-19 


N s = 1.5493*10 


Telescope area 

^ ^ tel “ ^obs \ 

re- =0.5871 

* / 


.15 


E r - 10 = 2.9003 


N g*T) -■ 5.4226* 10 


PROBABILITY OF MISSES AND FALSE ALARM 
The means of the APD outputs: 

noise only P 0 := ( n b + - — t] p 0 =53.5488 


1 . 6- 10 


signal and noise p i (n := n s + p . q 

The variance of the preamplifier noise and APD surface leakage current noise: 


var := 


N -t 
1 amp x 


I s T 


(i.eio 19 ) 1610 ' 


.-19 


:= ^var 


— = 1.7895- 10 2 
G 2 


Noise contributions 
from the preamp 
noise and APD 
surface leak current 


N -x 
amp t 


1 


' - 19 \^ G 

, 1.6 10 j 


= 13.3747 


l s x 1 


1 . 6 - 10" 19 G 


= 0.252 
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The p.d.f. of the APD output using Webb's approximation: 
APD excess noise factor 

F:=k eff G^ 2 -ij.(l-k eff ) F = 3.4734 

signal absent, noise only 


S 00 := J G F H 0 

PpDo( z ) := 


s go = 2 . 0457- 10 


1 

1 



-z 


G-(F-l)zj 

1 -t 

' VAL/ 

3 

2 

2* 

", G(F- l) z 

1 -I" 

s 00 


'00 


signal plus noise 


s ll( n s) : =J G F Hi(n s ) 


P PDl( z,n s) := ~F=" 

J 2 k 


1 + 


G (F - l)*z 


•exp 


- 1 

2 - 

G(F 

- l)-zl 

S 11 

S) 

II 


11 


(“•) 


Probabilities of false alarm: 

Calculate the Integration upper and lower limits 

-Gp 0 


lower limit of the integral 
Upper limit of the integral 


• 0 : 


•00 


z 1 := 100 


Standard deviation of the total noise 


s o := J s oo 2 + (j2 


s 0 = 2 . 8655-10 
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Solution for the normalized threshold level at a given false alarm rate 

M :* 10 k := 0, 1 .. M TOL := 10 12 The integration error tolerance 

NT k := 6.5 +■ 0.25 k Normalized threshold level wrt to the standard deviation of the total noise 

P fa k := P fa ( NT ») p FAk := i _ exp|~-pfa k J 


Find the point at PFA=0.01 
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The link margin calculation: 

N := 10 i := 0..N TOL := 10 5 

ns := 0.015 + !■■!— normalized signal level 

1 N 


PDj := P D( ns i' N s' 1 ! APD’^opt) 


F=yg 




- .98 

=#= 





f 





— I 






r 





_ 

f 
























\ 













0 0.05 0.1 0.15 0.2 0.25 

ns. 

t 


‘ 0.015 


6.234 1 • 10 3 

0.0175 


0.0237 

0.02 


0.069 

0.0225 


0.1597 

0.025 


0.3012 

0.0275 

PD = 

0.477 

0.03 


0.6535 

0.0325 


0.7988 

0.035 


0.898 

0.0375 


0.9549 

0.04 


0.9826 


Signal level at PD*(1-PFA)<=0.98 

PPD(x) := linterp(ns,PD,x) form a function of PD vs. normalized threshold 

x g := 0.013 initial guess 

sig ^ := root^PPD(x g ) - g j 

Link margin: 

sig ^ = 0.0407 10-log (— — J = 13.9085 (dB) 

\ s *8 min/ 


Link margin vs. transmitted laser pulse energy 

1 := 1 .. 20 

sig, := 0.05T SIG, := sig, E t 

I Sig l 1 

LM1. := 10- log 

\ si S min/ 
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0 10 20 30 40 50 60 70 80 90 100 


SIGj-10 3 

Laser Xmit Pulse Energy (mJ) 

Link margin vs. received pulse width: 
k := 1,2.. 80 
TFWHM k := k- 20- 10' 9 

Actual received signal vs. received pulse width 
(Hare the receiver was approximated as an integrator, 
since the distribution of a filtered APD output i9 too 
difficult to model at this time.) 

smax, := iff TFWHM. <x, 1 , } 

1 k TFWHMj 

/ smax, 

LM2 k := 10 log - 

min^ 


LM2 k 
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Link Margin vs. One-way Atmosphere Transmission 

minimun detectable signal if transmissions 00%: sigO ^ := sig ^Yri atmo 


M := 30 m := 1 .. M 


TATMO := — 
m M 


LM3 m := 10-log 


(TATMOj 2 

min 



o 0.5 1 

TATMO 

m 

WRITEPRN ( LMvsEt ) := augment(SIG.LMl) 

WRITEPRN ( LM vsFWHM ) := augment (TFWHM.LM2) 
WRITEPRN ( LM vsTATMO ) := augment ( TATMO, LM3) 
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Notes to the Mathcad Program 
for Laser Altimeter Link Analysis 

Xiaoli Sun, Code 924/JHU, May 18, 1995 

1. False Alarm Probability 

The average rate of false alarms can be calculated as 1 


where pf a is the false alarm probability at a fixed time and xw is the 
average noise pulse width above the threshold. It may be assumed that 
tw=1/BW3cjb- The number of false alarms over a given range gate time, 
Tq, follows a Poisson distribution and the false alarm probability is 
given by 


PpA 1 ex p| Tp A j = 1 ex P[— ToBWsjjBPfa] 

The false alarm probability at a fixed time can be written as 


(2) 


Pfa=| p(yin 0 ) d y= 

J y T 


>/2k<5 




PD^Wdxdy 


/y T 


(3) 


where p(yljxo) is the p.d.f. of the APD preamplifier output and is equal 
to the convolution of the p.d.f. of the APD output and the p.d.f. of the 
preamplifier noise. The latter is assumed to be additive and Gaussian 
with mean, |i g , and standard deviation, a. Usually, the preamplifier 
itself has zero DC offset, (j. g =l s x/q with l s the APD surface leakage 
current, x the equivalent integration time of the lowpass filter, and q 
the electron charge. The variance is given by 

<j2= (N amp +2I s )BW n x 2 /q 2 (4) 

with Namp the one sided input noise spectral density of the 
preamplifier in (A 2 /Hz) and BW n the noise bandwidth of the lowpass 
filter. 


The p.d.f. of the APD output in number of photoelectrons is given by 2 


1 M. I. Skolnik, Introduction to Radar Systems, McGraw-Hill, New York, 1962, p. 31. 

2 P. P. Webb, R. J. McIntyre, and J. Conradi, "Property of avalanche photodiodes,” RCA 
Review, Vol. 35, pp. 234-278, June 1974. 



Pp D (xl(J. 0 )- 


l 


\J 2kG 2 F(I 0 


1 + 


(x-Gp 0 )(F-l) 


GFjx 0 


3/2 eX P 


/ 


- (x-Gp 0 ) J 


\ '»r 2 Fu J I i < x - G ^o)(F-i) 

l 2G GF^ 


(5) 


where |io is the average number of primary photoelectrons given by 


APD^b 1 . fb X 

fi 0 - hf~ + q 


(6) 


with tiapd the APD quantum efficiency, Pb the received background 
noise power, hf the photon energy, and lb the APD bulk leakage current. 


Interchange the order of integrations in Eq. (3) and substitute 
u=(x-y+|i g )/a, and z=(x-Gpo)/(G 2 F|io) 1/2 . 


Pfa~ 


J2kg 


eX P 


(y-x-n,,) 2 ] 1 


yr 


j 


P PD( xl l I o) dx 


J y= ex p[- ^]du|p PD (xl^ 0 )dx 


/o 2 f ^ 0 


7^ exp t" J ^l du 


1 

1 

CXD ( -z 2 ) 

V2i 

L ( G(F-l)z 1 

3/2 ex Pi r G(F-l)zl 


v'G^Hol 

ITVo^Ji 


A normalized threshold is defined as 

yr^Gp-o+M-g) 


Y t = 


s/ G 2 F^ 0 +a 2 


(7) 

(8) 

(9) 


( 10 ) 


that is, Yj is how many standard deviations the actual threshold is 
above the average DC component. Usually, jig«Gjio, we can drop jig in 
the calculations. 


Substitute (10) into (7) 



| P PDz( zl Po) dz 


(11) 


Pfa= 


<9 


\j G 2 F)i 0 z-Y T \/ G 2 F)i 0 +a 2 


2. Probability of Correct Detection 

The probability of correction detection can be written as 


3 d=| P( x| Pi 

Jy t 


)dx= 


Yt 


f ^2ia exp [ K D(xlp i )dxdy 

*/ 0 


( 12 ) 


where |ii is the average number primary photoelectrons when the 
signal is present, 


Hi=Mo+n s (13) 

with n s the average number of detected signal photons per pulse. 

Similar to (7), interchange the order of integrations in Eq. (12) and 
substitute u=(x+^ g -y)/o, and z=(x-Gm)/(G 2 Fm) 1 /2 , 



where Yj is the normalized threshold defined in (10). 
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SL6140 

400MHz WIDEBAND AGC AMPLIFIER 

(Supersedes edition in September 1988 Linear 1C Handbook) 


The SL6140 is an integrated broadband AGC amplfier, 
designed on an advanced 3-micron all implanted bipolar 
process. The amplifier provides over 15dB of linear gain 
into 50Q at 400MHz. 

Accurate gain control is also provided with over 70dB of 
dynamic range. 

The SL6140 provides over 45dB of voltage gain with an 
R L of1kQ. 

FEATURES 

■ 400MHz Bandwidth (R L = 50Q) 

■ High Voltage Gain 45dB (R|_ = IkQ) 

■ 70dB Gain Control Range 

■ High Output Level at Low Gain 

■ Accurate Gain Control 

■ Full Military Temperature Range (CM only) 

■ MCI 590 Replacement with Improved Performance 

APPLICATIONS 

■ RF/IF Amplifier 

■ High Gain Mixers 

■ Video Amplifiers 


ORDERING INFORMATION 

SL6140 NA MP SL6140 B CM 

SL6140 A CM SL6140AC CM 
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ELECTRICAL CHARACTERISTICS 
Test conditions (unless otherwise stated): 

Tarnb""25*C, V(^= + 12V 


Characteristic 

Din 

Value 


Conditions 

“in 

Min 

Typ 

Max 

Units 

Supply current 

5.6,7 

♦ 

19 

23 

mA 


Output stage current 

5,6 

5 

7 

9 

mA 



(sum) 






Output current matching 
(magnitude of difference of output 
currents) 

5,6 


0.5 

* 

mA 


AGC range 

2 

60 

75 


dB 

See Fig. 4 & Note 1 

Voltage gain (single ended) 

5.6 

40 

45 

- 

dB 

Rl= Ikfl See Fig. 5 & Note 1 




55 


dB 

Tuned input and output 


5.6 

- 

15 

- 

dB 

R l = son 

Bandwidth (-3dB) 

5,6 

- 

25 

- 

MHz 

R(_ * 1 kO See Fig. 5 




400 


MHz 

R l = 500 

Maximum output level (single ended) 







OdB AGC 

5.6 

25 

3.5 

- 

V p-p 

Note 1 

-30dB AGC 

5.6 

2.5 

3.5 

- 

Vp-p 

R L = IkO, Note 1 


NOTE 1. Guaranteed but not tested for MP package 


DESCRIPTION 

The SL6140 (Fig.3) is a high gain amplifier with an 
AGO control capable of reducing the gain of the amplifier 
by over 70dB. The gain is adjustable by applying a 
voltage to the AGC input via an external resistor (Ragc)- 
the value of which adjusts the curve of gain reduction 
versus control voltage (see Fig. 4) As the output stage of 
the amplifier is an open collector the maximum voltage 
gain is determined by R L . With load resistance of 1 kO the 
single ended voltage gain is 45dB and with a load 
resistance of 50Q the voltage gain is l5dB (20log, 0 
Vout^in) Another parameter that depends on the load 
resistance is the bandwidth: 25MHz for R L = ikO, as 
compared with 400MHz for R L = 500. R L is chosen to give 
either the required bandwidth or voltage gain for the circuit. 

Fig. 7 shows the input impedance of the device. 
Accurate impedance matching to both inputs and outputs 
of this device (by resonant circuit or other means) can 
raise the gain to 55dB but for most circumstances a 500 
source impedance is adequate. 


ABSOLUTE MAXIMUM RATINGS 


Supply voltage, Vqq 

+ 18V 

Input voltage (differential) 

+ 5V 

AGC Supply 

Vqc 

Storage temperature 

-55°C to+ 150*C 

Operating temperature 

SL6140 MP 

O'C to 70*C 

SL6140 B CM 

-40*C to + 85*C 

SL6140 A CM 

-55 *C to+ 125*C 

Chip operating temperature 

SL6140 MP 

+ 175*C 

SL6140 (CM variants) 

+ 1 50 * C 

-IERMAL RESISTANCE 

Chip-to-ambient 

SL6140 MP 

225*C/W 

SL6140 (CM variants) 

163’C/W 

Chip-to-case 

SL6140 MP 

65° CAW 

SL6140 (CM variants) 

57*C/W 
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GLAS Master Oscillator Frequency 
Determinization 

Christopher T. Field 
October 17, 1996 


The GLAS instrument is a laser altimeter designed to measure the satel- 
lite altitude above the ground by measuring the time for a light pulse to 
travel to the ground and back. The range is to be measured to 7 cm from 
the 600 to 700 km orbit. There are many factors in the error budget; I allo- 
cate 0.7 cm to clock frequency error. Therefore, the frequency of the GLAS 
timing oscillator must be known to 7 mm/700 km or one part in 10 8 . The 
GLAS mission is expected to last 5 years so the oscillator frequency must 
change slower than one part in 5 x 10 8 per year or there must be a way to 
account for the change in frequency. This performance is near the limits of 
state of the art crystal oscillator performance. 

This report examines some of the factors involved in using on board GPS 
to provide time information which can be used to measure or control the 
master oscillator frequency. If sufficiently accurate frequency information 
can be determined, then the master oscillator may be of ordinary stability 
which decreases its cost and increases its reliability. 

In addition to knowing the master oscillator frequency, the position of 
the laser spot on the ground must be known to 3 meters. This requires that 
the position of the exit aperture must be known to a fraction of a meter at 
the moment the laser fires. The satellite is moving at about 7 km/sec, or 
7 m/msec. Assuming the orbit is known, knowledge of the time at which the 
laser fires to 0.1 msec enables precise knowledge of the exit aperture position. 
Consequently it would be helpful to synchronize (phase lock?) the master 
oscillator to GPS time. 

Tracking four GPS satellites permits time and position to be determined, 
called 3D mode. If the receive’s location is known, then only one satellite is 
required to determine time. This is known as ID mode. Additional satellites 
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over determine the solution and permit a least squares solution with a smaller 
position and/or time error variance since most errors from each GPS satellite 
are independent . Each satellite dithers its signal to provide time and location 
uncertainty. The satellite dithers are independent of each other and has a 
power spectrum extending from about 1 Hz to 0.5 mHz (1/2 hour period). 
Including extra satellites in the solution reduces the time and position effects 
of the dither. From the ground, between 4 and 7 satellites are visible. From 
low earth orbit, Tom Clark says at least 12 GPS satellites should be visible. 
While the dithers are not known in real time, their errors can be further 
reduced by post processing on the ground. 

Tom Clark has developed the Totally Accurate Clock (TAC) to replace 
or augment atomic clocks in remote sights. The TAC is a small GPS re- 
ceiver (Motorola’s ONCORE (PVT-6) OEM GPS receiver) which generates 
a pulse once each second. ONCORE is a six channel receiver which means it 
can track six satellites. Tom has done experiments comparing the difference 
between the GPS 1 PPS and the 1 PPS generated by various atomic stan- 
dards. It is known that the ONCORE has an inherent jitter of ±50 ns so 
Tom computes an average over 30 sec of the time difference between the PPS 
signals. In 3D mode, the 30 second average GPS 1 PPS has a peak jitter 
of ±500 ns around perfect 1 PPS timing. When in ID mode, the 30 second 
average difference drops to ±100 ns peak with an RMS error less than 40 ns. 

There are several ways to incorporate GPS time information for oscillator 
frequency measurement. The GPS 1 PPS may be used to measure either the 
master oscillator frequency or its period. Additionally the GPS signal may 
be used to steer the master oscillator to the desired frequency. Frequency 
measurement can be done by using the GPS 1 PPS to gate a running counter, 
one which is never reset, which counts cycles of the master oscillator. Each 
latched value is reported to the ground. The difference between consecu- 
tive readings is the number of oscillator cycles between the readings (in one 
second) from which the oscillator frequency can be determined. The aver- 
age frequency over a given, sufficiently long period is easily determined. To 
measure the period, the master oscillator can be divided to yield 1 PPS. The 
period can be determined from measurement of the time interval between the 
master PPS and GPS PPS. Active frequency control, provides a frequency 
feedback to change the oscillator frequency as it deviates from the desired 
rate. Ideally information on changes to the control signal are provided to the 
ground along with the satellite status information. 

The GLAS instrument requires several clock signals. There is an 800 MHz 
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signal to trigger the analog to digital converters, 50 MHz for the range gate, 
25 MHz for the digital signal processor, 5 or 2.5 MHz for the cloud channel, 
and 40 Hz for the laser trigger. These signals will be derived synchronously 
from one master oscillator by a digital divider chain. The general plan is to 
use the GPS 1 PPS and the divider chain to measure either the oscillator 
frequency or period. 

Before examining particular ways to measure frequency, here are some 
general definitions. Let the master oscillator have frequency f a . The GPS 
1 PPS has time jitter which makes the rising edge of the 1 PPS occur at 
times t Q , ti, • • ■. Ideally, the times are given by U = iAT where AT = 1 sec., 
but in fact are given by U = iAT 4- <5j where 6* is the jitter of the z’th edge. 
It will be assumed that <5 t < AT, Gaussian distributed with zero mean and 
statistically independent. 

For the frequency measurement approach, let the master oscillator drive 
a running counter, one which is never reset. At the leading edge of the 
GPS PPS, the counter value is latched to yield a value n. Let the set of 
observations no, nj , • • • , tim be made at times to, t\, • • • , t m- Assume that the 
master oscillator frequency, / 0 , is constant but unknown. Then 

n* = n 0 + (iAT + Si — 6 0 )f o (1) 

where <5j is zero mean Gaussian distributed with variance erf. Then the 
difference between the counter values at observation i and j is given by 

n t -n j = [(i-j) + S i -6 j ]ATf 0 (2) 

The estimated frequency during a single interval, (U,tj) is given by 

= 


By inspection, the expected value of is given by E = fo and 

(3b) provides an unbiased [1] estimate of the master oscillator frequency 
during the time interval. The variance of the estimated frequency is 


°l = {/•(«)} 

(4a) 

2cr J f 2 

(i- j) 2 AT 2 Jo 

(4b) 


7^1 ■j 


(i - j)AT 

6i - 8j 


1 + 


(i - J) AT, 


fo 


(3a) 

(3b) 
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which shows that the frequency estimate standard deviation drops as the the 
observation time interval increases. 

Counter values will be recorded each second for M seconds which yields 
M' = M+l counter readings which can be processed in several way. Adjacent 
counter values can be used to compute M one second frequency estimates 
which are averaged or (M + l)/2 frequency estimates can be computed from 
pairs and those values averaged. As shown below, for reasonable value of M 
the second approach yields a lower variance for the frequency estimate. Let 
the AT counter readings be made at times 6 0 , 1 + 6 X , • • • , M + S M . Between 
adjacent readings, a frequency estimate is made. The average estimate of 
frequency over the time t 0 to t M is given by 


1 y- rij - n,_x 

M “ AT 


(5 a) 


But because the counters are not reset after each reading all of the n values 
cancel except the first and last so 


fo 


r Hd_ — np 

MAT 



6m — 6q \ 
MAT ) 


fo 


(5b) 

(5c) 


As before, because the 8 s are zero mean this is an unbiased frequency esti- 
mate. The variance of the frequency estimate is computed as follows. 



1 yy Tlj - Wj-i 

M “ AT 



(6a) 


Again, all of the n’s except the first and last cancel so 


2 r 9 

a fo - (MMy fo 

Only the first and last count are needed to perform the average and the 
intermediate readings don’t reduce the estimate variance. 

An alternative approach permits contribution from all of the count read- 
ings and permits a smaller variance for the frequency estimate. Suppose that 
M' = M + 1 is even. The readings can be grouped into M' f 2 pairs spaced 
Ax AT apart. Then there are AT/2 estimates of frequency each given by 
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(rci'+Ai —J v)/ ( Ai AT) where the t's do not have to be time ordered. For exam- 
ple for A* = 1, the i' = 1, 2, 3, ■ • • , M'/2 correspond to i = 0, 2, 4, • • • , M - 1. 
The average estimated frequency is given by 


f 1 y? n (i / + Ai) - ?v 

io - M'j2 AiAT 


t'=i 

M'/2 


M'j 2 


E 1 + 


<5(»'+At) ~ 

AiAT 


/ 0 


(7a) 

(7b) 


Once again, (7b) is an unbiased estimate of the frequency. However, because 
each observation is included in the sum only once, the estimate variance, 
may be smaller than before. Taking into account that the Si s are zero mean 
and statistically independent, the variance is given by 



I [ 1 ^ W«>+At) ~ 

\ M'j 2 AiAT 

1 l f 2 

M'/2(AiAT) 2Jo 


2 ) 



(8a) 

(8b) 


The interval A i must be between 1 and M’j 2. Clearly using the largest 
possible value of A i produces the lowest variance given by 



8 2 a\ 2 

M' (M'AT) 2} ° 


(8c) 


Therefore for M' > 5 using all of the readings to produce M'j 2 frequency 
estimates provides a lower variance and hence a more efficient estimate of 
the oscillator frequency than using the first and last readings only. 

An alternative approach measures the master oscillator period. Assume 
f a /N is nearly 1 Hz. The terminal count from a divide by N counter will 
also be nearly 1 PPS. Call this signal the master PPS. See Figure 1. At time 
ti = iAT + Si , let AU + <5 At* be the measured time between the master PPS 
and the GPS PPS where A U is the true time interval and S Au is the error in 
the interval measurement. The period of the master PPS is P a = Nj f 0 so, 
at time the true time interval will be A U = (U - t 0 ) — k(i,o)P 0 /N and the 
measured time interval will be At* 4- 8& ti = (U — to) + ^At< — k(i, 0 )P o /N where 
k m is the number of master oscillator cycles between t t and t 0 . Generally 
= i but it may be different because of cycle slips. An estimate of the 
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Figure 1: Picture showing the time relationships between 1 PPS from the 
GPS and the master oscillator divide by N counter. 


average master oscillator period is given by 

p - ~ *») + (Atj — A ti) 

■* < r N 


(9a) 


_ j)Ar + (gj £j) + (Atj - Afj) + (6 At} — <^Atj ) jy 

^0.0 

where = Ar^o) — %,o)- Because the <5’s are zero mean and assuming 
that 6 Ati are also zero mean, this is an unbiased estimate of the period. The 
variance of the period measurement is given by 

. /\r2 

a p 0 = l 2cr s + 2cr At]r2— (10) 

K U.i) 

where a\ t is the variance of 6& t . 

As before, multiple time interval readings may be averaged in several 
ways. The straight average of P o(M _ 1} for * = 1, 2, 3, • • • , M yields the average 
and variance of the single pair of observations at time t 0 and t M as before 
and will not be investigated further. Using the alternative approach from 
before, the interval measurements may be paired with each value used only 
once. The average period is then be given by 

1 M '/ 2 _ 

P ° = W/2 (11) 

The variance of the estimate in (11) is calculated as before. 
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(12a) 


A = * {}*} - E 2 {L} 

— 1 2 aj + 2a% t jy2 

M' / 2 ^(j/ + At,t) 


(12b) 


If the frequency of the two oscillators are well matched and the timing jitter 
is not too large, fcq'+AM) can be bounded by 1 < %'+A«,i) < M'/2. Once 
again, the smallest variance is achieved by the largest value of k ~ M’ / 2. 
Therefore the variance of the period estimate is given by 



8 2g|+_2g|t 2 

M' M a ' 


(12c) 


The time interval between the two 1 PPS readings does not have to be 
measured with high accuracy because the required precision is provided by 
having a large value of i — j rather than high precision for At. The time 
interval is measured over a maximum of one second. Therefore, clock errors 
do not accumulate. If the timing clock is running fast by one part in 10 7 , the 
maximum interval error will be 100 nsec, less than the likely GPS PPS jitter 
in 3D mode. This error is likely to be similar for the interval at each end of 
the sampling time and cancel out. But even if it does not cancel, increasing 
A i will decrease the error in the master oscillator period measurement ar- 
bitrarily. Consequently AT* may be measured with a “low accuracy” clock 
such as the master oscillator itself and a separate oscillator is not required. 

The time interval is measured by starting a high frequency counter when 
the divide by N counter reaches zero and stopping (latching) it on the GPS 
PPS leading edge. But if the divide by N counter is a count up counter, it 
is zero when the master oscillator PPS is generated and its value when the 
GPS leading edge comes is the number of master oscillator cycles in the time 
interval. Therefore let rii be the value of the divide by N counter when the 
GPS leading edge comes. Then At* = (n* ± l)//o where the db accounts for 
the single count uncertainty. Using a 10 MHz counter, AT can be measured 
to one part in 10 7 and master oscillator drift of a few parts in 10 8 will not be 
important. Consequently no additional oscillators or electronics are needed 
to measure the lag time. This system then performs just the same as the 
previous one which latches the value of running counters. This system has the 
additional advantage that it always reports the time difference between the 
master oscillator, hence the laser firing time, and GPS time which simplifies 
determinization of the laser spatial position when it fired. 
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Glenn Lightsey and I have spoken about four space qualified GPS re- 
ceivers. The Jet Propulsion Lab (JPL) is developing the Turbo Rogu, Trim- 
bol Navigation has the Vector, Space Systems Loral has the Tensor, and 
Trimbol and Honeywell are developing a new, unnamed, receiver. 

The Turbo Rogu is a dual frequency receiver developed to probe the at- 
mosphere by measuring the ionic delay as GPS satellites move above or below 
the horizon. The Turbo Rogu does not have attitude measuring capability 
and was not intended as the principle navigation receiver. The Turo Rogu 
has flown on a few small missions. 

The Vector by Trimbol Navigation is a commercial development which 
uses commercial parts. It is a six channel receiver capable of reporting time, 
position, and attitude. It is a low cost unit, $30,000 and has flown 6 times. 
The typical orbits have been 800 km and lifetime was 6 to 18 months. Al- 
though the functions required for GLAS are there, the lifetime is not. 

Space System Loral took the Trimbol Navigation design and engineered it 
for space flight . The Tensor (do you see a pattern here? ) is a fully redundant 
9 channel unit designed to provide time, position, and attitude for 5 years 
through 100 krad radiation. The Tensor costs about $500,000. 

Trimbol in cooperation with Honeywell and GSFC Code 712 is developing 
another GPS receiver to cost about $100,000 and provide the functions of 
the Tensor but may not have quite the lifetime. 

Many GPS receivers, including the ones above, provide a 1 PPS output 
which can be used by the circuits described in this paper to measure a crystal 
oscillator frequency. In orbit, the rising edge of the Tensor’s 1 PPS signal is 
expected to be within l^sec of the absolute time point. With post processing, 
it is possible to reduce the uncertainty of the edge location. With 1/isec 
timing uncertainty, the observation time must be 100 to 1000 seconds long 
to provide frequency or period knowledge to one part in 10 8 . 

I thank Bryan Blair, Tom Clark, and Glenn Lightsey all of Goddard 
Space Flight Center for providing extensive background information on GPS 
receivers and position determination. 

More information on the TAC can be found on the Wide World Web at 
ftp://aleph.gsfc.nasa.gov/GPS/totally.accurate.clock/ 
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